Manipulating dipolar and spin-exchange interactions 
in spin-1 Bose-Einstein condensates 
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It remains a challenge to independently manipulate the magnetic dipolar and the spin-exchange 
interactions, which are entangled in many spin systems, particularly in spin-1 Bose-Einstein conden- 
sates. For this purpose, we put forward a sequence of rf pulses and the periodic dynamical decoupling 
sequence of optical Feshbach resonance pulses to control the dipolar and the spin-exchange interac- 
tions, respectively. Our analytic results and the numerical simulations demonstrate that either of 
the two interactions can be suppressed to make the other dominate the spin dynamics; furthermore, 
both of the interactions can be simultaneously suppressed to realize spinor-condensate-based mag- 
netometers with a higher sensitivity. This manipulation method may find its wide applications in 
magnetic resonance and spintronics. 

PACS numbers: 03.75.Mn, 03.75.Kk, 03.67.Pp 



Introduction — Long-range magnetic dipolar interac- 
tion and short-range spin-exchange interaction are fun- 
damental in many spin systems, such as the electron 
and nuclear spin systems [l| , the ferromagnetic ultrathin 
films Q , the ultracold atomic and molecular gases [!, 0] , 
and so on. The competition between the two spin in- 
teractions diversifies not only the ground state but also 
the spin dynamics, especially in spinor &11| or dipo- 
lar Bose-Einstein condensates (BECs) jl2|-|l8l|. where the 
internal spin deg rees of freedom are released in opti- 
cal traps fl9L Ho] and many physical parameters can 
be tuned precisely. However, the entanglement of the 
dipolar and spin-exchange interactions makes it difficult 
to understand the observed phenomena in many exper- 
iments. Such an interesting example is the spin texture 
or the spin domain in 87 Rb spin-1 condensates, which 
might be ascribed to either one of these two interactions 
or both of them [U [22[. Accordingly, it is highly de- 
sired to suppress one spin interaction and make the other 
dominant. Although efforts have been made to prohibit 
either the dipolar interaction [2l[ or the spin exchange 
interaction 0, [24| , it is yet difficult to clearly separate 
the individual effects. 

In this paper, we propose to independently suppress 
the dipolar and the spin-exchange interaction by a se- 
quence of rf pulses ( Fig. [T]) and optical Feshbach pe- 
riodic dynamical decoupling (PDD) sequences (24j |. re- 
spectively. This scheme was verified by the analytic 
derivations, demonstrating that the manipulations of the 
relative strength between the two spin interactions en- 
able experimentalists to unambiguously distinguish the 
specific role of the dipolar and the spin-exchange in- 
teractions. It is notable that the sequence of rf pulses 
is similar to the so called WAHUHA rf pulse sequence 



used to remove the secular dipolar interaction in nu- 
clear spin systems [25|, but our analytic results show 
that the WAHUHA sequence fails in spin-1 BECs while 
the pulse sequence shown in Fig.l, named as generalized 
WAHUHA (g- WAHUHA) rf pulse sequence in the fol- 
lowing text, works well in any spin dipolar system. To 
verify the manipulation means under specific experimen- 
tal conditions, we carried out numerical simulations with 
parameters of 87 Rb spin-1 condensates and the results 
are in agreement with our analytic conclusions. 

The model — Within the mean field theory, coupled 
Gross-Pitaevskii equations for a spin-1 BEC in an arbi- 
trary trap are[ 
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where M is the atom mass, ip a (a = ±1, 0) are the order 
parameters for the three components, N = i^a^a 
is the total number density, S — ^ Q ^ ip^F a ^ipp is 
the spin density with F the atom spin matrix, and the 
dipole integral operator is -Dr,(r) = J dr' r^jrs [^( r/ ) — 
3e^S(r') • e] (77 = x, y, z) with e the unit dipolar vector; 
S± = {S x ± iS y )/V2 and D± = {D x ± iD y )/V2. The 
spin-independent interaction and the spin-exchange in- 
teraction coefficients are, respectively, cq — 4irh 2 (ao + 
2a2)/3M and c 2 — Anh 2 (a 2 — ao)/3M with 00,2 denot- 
ing the s-wave scattering length in the two symmet- 
ric channels; the dipolar interaction coefficient is c d = 
MoSfMb/^ with no the vacuum magnetic permeability, 
fj,B the Bohr magneton, and gF the Lande (/-factor. 

We proceed by adopting the single mode approxima- 
tion (SMA) [28-30] . which separates the spinor conden- 
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FIG. 1. Generalized WAHUHA pulse sequence of one period. 
The vertical lines mark the time of each applied rf pulse and 
the alphabets above the vertical lines represent the axis along 
which the spins are rotated tt/2. 



sate wave function into a spin-independent spatial part 
<j>(r) and the spin part £, i.e., ^> a (r) — <fi(r)£ a . By drop- 
ping the density dependent terms, which do not affect the 
spin dynamics, we obtain the Hamiltonian under SMA 

Usma = c' 2 f + c' d (3f z - f) + 3cgC# - /J), (2) 

where the isotropic spin-exchange term is 
c 2 = ( c 2/2) / drp 2 (v) with p(r) = |0(r)| 2 , the anisotropic 
dipolar term contains two parts, the cylindrically sym- 
metric one c' d = (cy/4) / drdr' \ r J r , \ i p(r)p(r')(l - 

3cos 2 # e ) and the nonsymmetric one c d = 
(q/4)/ dvdv' ¥ J FW p(r)p(r') sin 2 9 e e ±i2 ^ with 9 e 
and (fie being the polar and azimuthal angles of (r — r'). 
The amplitudes of c' d and c' d are approximately in the 
same order but are one order of magnitude smaller than 
that of c 2 for 87 Rb condensates [31(. The spins are 
f x ,y, z = (t;\F x , v , z \0 and p = ft + / 2 + fl The spin 
dynamics under SMA becomes 

d 

±A 2 (n +1 - n_i)£ ±1 

+M$(C+i + Ci) + + C-i)], (3) 

+A 3 [Coi€+i +€-ii 2 + +c-i) 2 ], 

where A 1 = 2(c' 2 -c' d + 3c' d '), A 2 = 6(c' d - c£), and A 3 = 
— 6c^'. We have defined £ a = y/n^e lSa: with n a and <5 a 
being the fractional population and the phase of the crth 
component, respectively [Hj]. 

Independent suppression of the dipolar interaction — 
The Hamiltonian in Eq.([2]) reminds us two widely 
used methods, the magic angle and the WAHUHA se- 
quence [lU , which are capable of suppressing the dipolar 
interaction in nuclear spin systems. However, both of the 
two methods fail in a spin-1 BEC in that the nonsym- 
metric term with c d in Eq.© cannot be canceled [34|. 

We then generalize the WAHUHA sequence by adding 
extra 7r/2-rf pulses along x direction (see Fig. [T]). Each 
pulse instantaneously rotate the spin 7r/2 along the des- 
ignated direction. Due to its isotropic form, the spin- 
exchange term is not affected by any spin rotation, which 
provides us the opportunity to suppress the dipolar in- 
teraction independently. Below we show that this g- 
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FIG. 2. (Color Online) Time evolutions of (a) no and (b) 
fz- The black dashed lines denote the free evolution with 
both spin-exchange and dipolar interactions, the blue solid 
lines denote the free evolution without dipolar interaction, 
and the red crosses denote the evolution under a sequence 
of 100 g- WAHUHA periods. The coincidence of the crosses 
and the solid lines clearly indicates that the dipolar interac- 
tion is suppressed by the applied g- WAHUHA sequence with- 
out affecting the spin-exchange effect. The initial values of 
(n+i,no,n_i) are (0.4, 0.5, 0.1) and the geometry of the trap 
is chosen as c' d = 0.2c 2 and c d = — O.IC2 with c' 2 = — 1. 

WAHUHA sequence indeed prohibits dipolar effects for 
an arbitrary shape condensate. 

We define the free evolution during interval t £ 
[tit, (n + 1)t] between two adjacent pulses as Uq = 
exp(—iHsMAT~) and the rotation along the designated 77 
axis as U n = cxp(—iF v ~). The evolution for a complete 
g- WAHUHA period can be split into two halves, Ut = 
U-U+, where U± = U U x U U- y U U± x U U y U U- x U& 
In the toggling frame which is defined by the rf pulses [25], 
[Hj], the first half-period evolution operator becomes 
U+ = U 6 U 5 UaU 3 U 2 Ui where Uj = exp(-ifl J -r) (j = 
1,2,- •• ,6) with 
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According to the average Hamiltonian theory [36| . 
the evolution operator U+ can be reexpressed as U+ = 
exp(-ifff) with H = T,kLo H(k) - The term H(k) is 
proportional to r k and can be obtained with Magnus ex- 
pansion. In the limit of short r, the lowest order term 
H( > , which is independent of r, becomes dominant 

ff (0) = ^I> = 4/ 2 - (5) 

i—l 

Higher order terms are neglected since r is small. In 
current experiments [2l| . where the condensate density 
is about 10 14 cm~ 3 , c' d ~ c" d ~ 0Ac' 2 ~ 1Hz, the pulse 
delay is r ~ 10 _3 s, so that the higher order terms can be 
safely neglected, - 1Q- 4 H^\ 
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FIG. 3. (Color Online) Same as Fig. [2] except that the blue 
solid lines denote the free evolution without spin-exchange 
interaction and the red crosses denote the evolution under 
a sequence of 300 PDD pulses. The coincidence of the blue 
solid lines and the red crosses in both (a) and (b) indicates 
that the spin-exchange interaction are suppressed by the PDD 
sequence without affecting the dipolar effect. In (b), the black 
dashed line coincides with the blue solid one, manifesting the 
fact that the spin-exchange interaction does not affect the f z 
dynamics. 

The average Hamiltonian given by Eq. (JSJ clearly indi- 
cates that half of the g-WAHUHA sequence completely 
removes the dipolar interaction in the condensate to the 
lowest order r. Note that this result is universal, though 
the parameters c' d and c d closely depend on the conden- 
sate shape. Finally, the toggling frame does not coin- 
cide with the laboratory frame because of the extra rota- 
tion term U x U~ y U x U y U- x , introduced by the applied five 
pulses. By simply adding five more pulses (see Fig. [T]), 
the second-half rotation term counteracts the first-half 
one. In this way, the toggling frame exactly coincides 
with the laboratory frame after a complete g-WAHUHA 
sequence. 

To justify the above analysis of the g-WAHUHA 
method, we simulate the time evolution of a spin-1 sys- 
tem based on Eq.© in three cases: (i) free evolution 
with both dipolar and spin-exchange interaction; (ii) free 
evolution with only the spin-exchange interaction; (iii) 
evolution with both spin interactions and g-WAHUHA 
sequence. If the results for case (ii) and (iii) are close, 
we safely conclude that the g-WAHUHA sequence effec- 
tively and independently suppresses the dipolar interac- 
tion. We show in Fig. [2] no and f z , which are easily 
measured in experiments. As can be clearly seen, the 
red crosses coincide with the blue solid lines, indicating 
the perfect agreement of our previous analysis and the 
numerical simulations, i.e., the g-WAHUHA sequence in- 
dependently suppresses the dipolar interaction. 

Independent suppression of the spin- exchange 
interaction — We next consider the independent 
prohibition of the spin-exchange interaction in spin-1 
BECs, with the presence of the dipolar interaction. 
Previous works demonstrate that PDD pulse sequences 
are able to suppress the spin dynamics induced solely 
by c' 2 term in Eq.®, by flipping the sign of ci through 
the optical Feshbach resonance technique 0, [2H [37| . 
With the inclusion of dipolar terms, whether the PDD 



FIG. 4. (Color Online) Time evolution of (a) no and (b)/ z . 
The dashed lines denote the free evolution and the red crosses 
are under a sequence of 100 g-WAHUHA and PDD periods 
simultaneously. The freezing of the spin dynamics (crosses) 
indicates that both the dipolar and the spin-exchange inter- 
actions are suppressed. Parameters are the same as Fig. [2] 

sequences perform the same as above is unknown. We 
show in the following within the average Hamiltonian 
theory that the optical PDD sequences indeed suppress 
the spin-exchange effects, independent of the inclusion 
of the dipolar interaction. 

For a 2r PDD sequence, the evolution operator is Ut = 
U+U- where U± = exp(-iff ± r) with H± = ±c' 2 f 2 + 
c'AZfz - P) + Kify ~ /I); With the help of Magnus 
ex pan sion, the average Hamiltonian to the lowest order 
is [36| 

H^=c> d (3f z -f) + 3ci>(f y -f x ), (6) 

which is free of c' 2 and r, and includes only the dipo- 
lar terms. Therefore, in the limit of small r, the op- 
tical PDD sequence removes the spin-exchange effects 
and leaves only the dipolar terms in a spin-1 conden- 
sate. Note that with a pure spin-exchange interaction, 
the PDD sequence can remove all the average Hamilto- 
nian terms [2j]. Numerical results from Eq. © for no 
and f z are shown in Fig. [51 As one can see, the time 
evolution without spin exchange interaction (blue solid 
lines) is the same as that under the PDD sequence (red 
crosses), which confirms evidently the suppression of the 
spin-exchange interaction by the PDD sequence without 
affecting the dipolar effect. 

Simultaneous suppression of both dipolar and spin- 
exchange interactions — In some applications, for in- 
stance, a magnetometer based on spin-1 BECs [38|, it 
is desired to simultaneously suppress both spin interac- 
tions: the spin-exchange and the dipolar one. Since the 
g-WAHUHA sequence with rf pulses and the PDD se- 
quence with optical ones are independent, it is straight- 
forward to make a cooperative sequence which applies 
the rf pulses and the optical ones at the same time. For 
convenience, we set the period of the g-WAHUHA se- 
quence to be the same as that of the PDD sequence. We 
present in Fig. [4] the frozen spin dynamics of the spin- 
1 condensates where the dipolar and the spin-exchange 
interactions are both suppressed. 

Beyond the single mode approximation — Our analysis 
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FIG. 5. (Color Online) Numerical simulations beyond SMA 
for a spin-1 BEC with spin-exchange and dipolar interac- 
tions. The system is either free or under the g-WAHUHA 
sequences. In (a): Time dependence of the fraction of each 
component, n+i (solid line), no (dashed line), and n_i (dot- 
ted line) in free evolution; n+i (crosses), no (circles), and 
n_i (asterisks) under the g-WAHUHA sequences [the pulse 
delay is r = 0.005(1/cj z )]. The freely evolved density distri- 
butions in the x — z plane at t = 30 are shown in (b), (c), 
and (d) for +1,0,-1 components, respectively. The density 
distributions under g-WAHUHA sequences are shown cor- 
respondingly in (e), (f), and (g). The SMA is invalid for 
a freely evolved spin-1 condensate but it is kept if the g- 
WAHUHA sequences are applied, indicating the suppression 
of dipolar interactions. The initial state is a ground state of a 
87 Rb spin-1 condensate with pure spin-exchange interaction. 
The total atom number is 2 x 10 5 , the initial fractions are 
n+i = 0.767, no = 0.218, n_i = 0.015, and the aspect ratio is 
uj x : u) y : oj z — 2 : 4 : 1 with uu z — (2irx) 20Hz. The dipolar 
coefficient Cd is amplified five times to violate substantially 
the SMA during the free evolution. 



the dipolar interaction. 

The simulation starts from a SMA ground state with 
pure spin-exchange interaction, where the spin dynamics 
can be only introduced by the dipolar interactions. We 
turn on the dipolar interaction at t = and deliberately 
amplify its strength by five times to test the power of g- 
WAHUHA sequences. Two numerical simulations based 
on Eq. |T]) are presented in Fig. [5] : In one simulation, we 
let the system evolve freely; in the other, we apply the 
g-WAHUHA sequences repeatedly. As shown clearly in 
Fig-EJa), the spin dynamics induced by the dipolar inter- 
action is evidently suppressed and the fractions of each 
spin component are frozen, though the SMA is signifi- 
cantly violated during the free evolution [see Fig. [5jb), 
(c), (d)]. 

Conclusion — The proposed g-WAHUAHA sequence 
of rf pulses and the PDD seuqnces of optical Feshbach 
resonance pulses provide an experimentally feasible way 
to separate clearly the contribution of the dipolar and 
the spin-exchange interactions, which are usually entan- 
gled together in spin-1 BECs. Furthermore, by cooper- 
atively suppressing both spin interactions, the precision 
of a magnetometer based on spin-1 Bose condensates can 
be significantly improved. Generally, the developed g- 
WAHUHA sequence for suppressing the dipolar interac- 
tion paves a new way to determine the atomic and molec- 
ular structures at zero or low magnetic fields in magnetic 
resonance imaging, nuclear magnetic resonance, chem- 
istry, and biology. This sequence may also be utilized to 
extend the quantum coherence time of electron or nuclear 
spin systems in spintronics and quantum computing. 
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the following that the g-WAHUHA sequences also sup- 
press the formation of spin spatial structure caused by 
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